Critical Behavior of the Random Potts Chain 
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We study the critical behavior of the random g-state Potts quantum chain by density matrix 
renormahzation techniques. Critical exponents are calculated by scaling analysis of finite lattice 
data of short chains (L < 16) averaging over all possible realizations of disorder configurations 
chosen according to a binary distribution. Our numerical results show that the critical properties of 
the model are independent of q in agreement with a renormahzation group analysis of Senthil and 
Majumdar (Phys. Rev. Lett. 76, 3001 (1996)). We show how an accurate analysis of moments of the 
distribution of magnetizations allows a precise determination of critical exponents, circumventing 
some problems related to binary disorder. Multiscaling properties of the model and dynamical 
correlation functions are also investigated. 
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I. INTRODUCTION 



The study of the critical properties of systems subject 
to quenched randomness is of great interest because of 
its relevance for many experimental systems. This is a 
challenging research area which has been quite active in 
the last decades. One has witnessed a growing interest 
in understanding the effect of randomness (bond and site 
disorder), since Harris showed that disorder is a relevant 
perturbation in a system undergoing a second order phase 
transition when the exponent a of the specific heat in the 
pure system is positive.EJ In this context the two dimen- 
sional Ising model has attracted considerable interest, 
since the Onsager singularity of the specific heat (a = 0) 
makes the Harris criterion inconclusive (for a.review, see 
e.g. Ref. y). Renormalizatioia, group studies, a supported 
by Monte Carlo simulations p established that homoge- 
neous disorder is marginally irrelevant, i.e. it does not 
modify the critical behavior of the system, except for the 
appearance of logarithmic corrections. 

In accordance with the Harris criterion, randomness is 
on the other hand a relevant perturbation for the q — ?i 
and g = 4 Potts models. Furthermore, an infinitesimal 
amount of disorder was shown to turn the first order 
phase transiti(wi|-pf the q > A Potts model into a sec- 
ond order one.EI^LJ Numerical estimates of the critical ex- 
ponents for an homogeneous disorder indeed showed a 
monotonic increase of the magnetization exponent (3 with 
the value q while the thermal exHcpent v remains con- 
stant within numerical accuracy.El"til Thus for the two 
dimensional classical Potts model, homogeneous disorder 
gives rise to several difi^erent universality classes depen- 
dent on the value of q. 

Other systems with quenched randomness which have 



been investigated in the past years are quantum spin 
chains. The interest in these models haS|-grown consid- 
erably after a remarkable paper of FishertS who derived 
many of the critical properties of the random transverse 
Ising chain (RTIC) using a real space renormalization- 
group scheme which is claimed to be asymptotically ex- 
act. His predictions were checked numerically by sev- 
eral authors, using the mapping of the RTIC onto a frea 
fermion modeljlj density-matrix rerpi'malization groupliil 
(DMRG) or quantum Monte CarlollJ Critical exponents 
calculated numerically were found in good agreement 
with Fisher's results. The renormahzation scheme of 
Fisher was also applied to the QpState disordered Potts 
chain by Senthil and Majumdar;t3 they found that the 
critical behavior does not depend on g, therefore all criti- 
cal exponents should be identical to those of the RTIC ob- 
tained by Fisher and corresponding to the case q — 2. We 
recall that the disordered quantum Potts chain is equiva- 
lent to a two dimensional classical Potts model with cop-, 
related disorder (analogous to the McCoy- Wu modelEII 
corresponding to q = 2), where the quantum formula- 
tion is obtained rfcom the classical one in the extreme 
anisotropic limit .t2l The conclusion of a unique universal- 
ity class is thus different from what found in the classical 
Potts model with homogeneous disorder discussed above. 
In this paper, we check the predictions of a unique 
universality class, independent of the number of states 
q, by studying bulk and surface magnetization proper- 
ties. We consider both the regimes g < 4 and g > 4, 
where the pure model exhibits a second and first-order 
transition, respectively. We also report an investigation 
of the multiscaling behavior of the magnetization, which 
strongly supports the asymptotic expression for the prob- 
ability distribution found by Fisher in the RTIC case. 
The dynamical behavior is considered through the de- 



* Unite Mixte de Recherche CNRS No 7556 



cay of the spin-spin correlation function along the time 
direction. All numerical results are obtaiaftdiby density 
matrix renormalization group techniques jESEa applied to 
the two dimensional classicaLPotts model with correlated 
disorder, ffie use Nishino'£l version of White's DMRG 
algorithm,llj where one renormalizes classical transfer 
matrices. The exponents are calculated from finite size 
scaling analysis of finite lattice data for chains of lengths 
L < 16 with exact enumeration over all possible (2^) 
disorder realizations, obtained from a binary distribu- 
tion. The justification for this type of approach will be 
presented in more detail in the text. 

Our numerical results are fully copsistent with the con- 
clusion of Scnthil and Majumdar,E£l i.e. the g-state Potts 
model has critical properties which do not depend on q. 
Moreover, dynamical properties which are not accessible 
by Fisher's renormalization group, do not seem to depend 
on q. 

The paper is organized as follows: In Sec. |l|, the model 
is presented and we give a justification for the choice of 
the numerical methods that are used. 



Ill 



is dc- 

Multi- 



Section 
voted to the estimates of the critical exponents 
fractality is investigated in Sec. IV. The dynamical spin- 
spin correlation functions are studied in Sec. |V|. Finally, 
a brief summary and an outlook will be presented in 



Sec. VI 



FIG. 1. Layered structure of the model with coupling con- 
stants Ki and K2 taken from the distribution (H) and indi- 
cated by thick and thin bonds, respectively. The arrows on 
the left indicate fixed boundary conditions, while spins on the 
right surface are left free. 



II. MODEL AND NUMERICAL METHODS 



We consider a randomly layered (see Fig. [|) two- 
dimensional classical g-state Potts model on a square lat- 
tice. The model is defined by the Hamiltonian 
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1^1,71,1^1 + 1,1 
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(1) 
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where the spins can take the values (T;^„ = 0, . . . , (7 — 1 
and the exchange couplings Ki are randomly chosen ac- 
cording to the binary distribution 



V{Ki^ 



1 



[6{Ki - Ki) + 6{Ki - K2)] . 



(2) 



Since the model is equivalent to the random xpantum 
Potts chain studied by Senthil and Majumdar,t3 we use 
the quantum chain's vocabulary. In the case q = 2, the 
system reduces to the RTIC. 

The particular arrangement of coupling constants, i.e. 
the fact that vertical couplings are identical to adjacent 
horizontal couplings to their right side (see Fig. m, and 
that both couplings Ki and K2 have equal probabilities 
(see Eq. (H)), makes the model self-j4ual. The location of 
the critical point is known exactly:E3 



„Ki 



1 



^K2 



-I =,. 



(3) 



The strength of the disorder is monitored by the ratio 
r = Ki/ K2, r — I corresponding to the pure model. 

Although the DMRG method is capable of dealing with 
very large systems (of few hundreds of sites) with remark- 
able accuracy, some care has to be taken when working 
with disordered systems. The infinite lattice version of 
the algorithm is known to give incorrect-results £3 es- 
pecially at strong disorder; it is essentiaE^ when deal- 
ing with inhomogeueous systems to use the finite lattice 
DMRG algorithmlij which is designed to determine ac- 
curately properties of finite systems (In our calculations 
we typically used three full sweeps through the lattice). 
It is furthewpore well known that due to the lack of self- 
averagingpJc3 an insufficient number of samples yields 
a wrong estimation of the average over randomness and 
might lead to twical values instead of the expected aver- 
age quantities .E3 As we aim to an accurate determination 
of critical exponents, we have chosen to restrict ourselves 
to relatively narrow strips (L < 16 which turns out to 
be large enough, since corrections to scaling will appear 
very small) for which we were able to do an exact enu- 
meration over all disorder realizations. As the system 
size is small compared to typical DMRG calculations a 
limited number of states kept {m < 16) is sufficient to 
obtain accurate results. 

The setup of the calculation follows Ref. ^^, where 
the DMRG method was used to calculate exponents and 
magnetization profiles for the pure Potts model. The 
Zq-symmetry of the Potts model has been broken by in- 
troducing fixed-free boundary conditions: We fixed the 
spin at one edge of the system to a value cro,n = 0, while 
at the other edge the spin is free. This type of boundaoz. 
conditions induces a finite magnetization in the systemEHl 
from which surface and bulk critical exponents can be 
calculated using finite size scaling analysis. 

Figure pi shows magnetization profiles (defined by 
Eq. (0)) of 10 samples for the three-state Potts model 



with L = 16 and with disorder amphtude r — 2. Esti- 
mated error bars are much smaller than symbol sizes; av- 
eraging over all disorder realizations produces the smooth 
profiles shown in Fig. y. In this procedure the only source 
of error comes from the truncation on the DMRG basis, 
this error is quite small even for to = 16 states kept as 
we restricted ourselves to small systems. 



s 




FIG. 2. Magnetization profiles for ten random configura- 
tions of exchange couplings at g = 3, r = 2, L = 16. Error 
bars are smaller than the symbols. 



III. CRITICAL EXPONENTS 

A. Finte-Size scaling 

We used the standard order parameter of the Potts 
model: 



NL(0]av = 



9[(^<T,,„,o)], 



9- 



1 



(4) 



where (...) denotes the thermal average and [. . .]^^ the 
average over randomness. Standard finite size techniques 
have been used to estimate the bulk and surface critical 
exponents 



'LJav 



L- 



'iJav 



where we have defined [to.^] 



LJav 



(5) 

[TOL(L/2)]av and 
[TTi^av = ['TiL(L)]av, which are the average magnetiza- 
tions in the middle of the strip and at the free edge. From 
finite size data we calculated the approximants, defined 
as 



Xb{L) 



ln[TO^ , Jav - ln[TO^ 



L-lJav 



ln(L + 1) - ln(L - 1) 



(6) 



in the bulk case and similarly for the surface approx- 
imants. In the limit L — > c» the approximants scale 
towards the bulk and surface exponents Xb and xi. 

Table | reports the finite size approximants for the sur- 
face exponents for various values of q and of disorder ra- 
tios r. The extrapolation towards the titermodynamic 
limit was done using the BST algorithm.123 The bound- 
ary magnetization exponent, xi, exhibits a striking con- 
vergence for all values of q and all disorder amplitudes, 
but for crossover effects discussed below. It is clearly in- 
dependent of q, and in good agreement with the RTIC 
value xi = 1/2. The best numerical results are obtained 
for g = 2, as the truncation error in the DMRG method is 
the smallest for this case, furthermore the surface expo- 
nent for the pure and disordered cases coincide for q = 2 
(they are both equal to 1/2), therefore no crossover ef- 
fects are to be expected. 



TABLE L Approximants of the surface exponent for various q and r. The last row shows the values extrapolated with the 
BST algorithm. We note that ai q — 6 the disorder amplitude being quite small, crossover effects from the pure model value 
a;i = 3 slow down the convergence. 





q 


= 2 




g = 3 




9 


= 4 


g = 6 


L 


r = 2 


r = 8 


r = 2 


r = 4 


r = 8 


r = 5 


r = 10 


r = 1.5 


5 

7 

9 

11 

13 

15 


0.41608413 
0.43773154 
0.45043349 
0.45880966 
0.46475540 
0.46919708 


0.44950692 
0.46307866 
0.47083067 
0.47587114 
0.47941816 
0.48205234 


0.46559328 
0.48454642 
0.49365387 
0.49851791 
0.50127872 


0.46393358 
0.47694537 
0.48353868 
0.48742293 
0.48993505 
0.49166836 


0.46844808 
0.47932515 
0.48475158 
0.48800690 
0.49014411 
0.49150054 


0.47671724 
0.48718127 
0.49161650 
0.49383119 
0.49481028 


0.47741488 
0.48819610 
0.49204019 
0.49387188 
0.49435771 


0.55719868 
0.56816316 
0.56481411 
0.56108433 
0.54974464 


oo 


0.50000(1) 


0.5000(1) 


0.498(4) 


0.499(1) 


0.500(1) 


0.496(10) 


0.50(5) 


- 



We recall that the values of the surface exponents for 
the pure Potts model are xi — 2/3 for q = 3, xi = 1 
ior q — 4 and xi — 3 for q > A. The last exponaat has 
also been recently calculated by DMRG methodS (the 



surface transition is continuous even if the bulk has a 
first order transition). In the case g = 6, r = 1.5, one 
notices that the approximants approach the value 1/2 in 
a non-monotonic fashion (the same situation occurs with 



g = 4, r = 2). This case corresponds to the weakest dis- 
order considered and the non-monotonicity is most hkely 
due to crossover effects from the exponent of the pure sys- 
tem (xi =3). It is not possible to extrapolate the data 
in the last column of Table | with the BST method, since 
the largest size analyzed is still in the crossover region. 
However we can give an estimate of the exponent using 
the analysis of the moments presented in Section IV B; 
from such an analysis we find 0.49(3) in good agreement 
with the value of the RTIC, 1/2. We also mention that 
calculations have been performed at g = 6, r = 4. The 
approximants then do not suffer from so strong crossover 
effects, an observation which supports the explanation 
given above. 

Table ^ shows the approximants of the bulk exponent 
obtained from Eq. (pf) for various q and r. As it can 
be seen the approximants show a non-monotonic behav- 
ior, with oscillations which are stronger as the disorder 
ratio increases. These-.oscillations, which were already 
observed in the RTIC,l3 are due to the choice of a binary 
distribution of exchange couplings which introduces an 
energy scale Ki — K2 in the problem; in other types of 
distributions_e.g. continuous ones, these oscillations are 
not present .E3 

Lattice extrapolation techniques, as the BST method, 
are not able to deal with oscillating corrections to scaling. 
We recall that the BST algorithm generates from an orig- 
inal series of iV approximants new series of A^ — 1, N ~2 



. . . approximants at each iteration which are expected to 
have faster convergence towards the asymptotic value. 
In the BST analysis we found that even when starting 
from a monotonic series of data at weak disorder ratios, 
as in the cases q = 2, r ^ 2, the next series generated 
by the algorithm do not converge at all, but they show 
oscillating behavior. Therefore the BST method turns 
out to be useless in the analysis of the bulk exponent 
data. A better way to get a precise estimate of bulk ex- 
ponent is based on the analysis of moments, which will 
be explained in Section IV E| . 

In any case the approximants of Table |^ for different 
values of q are rather close to each other, which supports 
again the idea of g-independence of their values. We re- 
call that for the pure Potts model the bulk exponent is 
Xb = 1/8 for g = 2, 4 and Xb = 2/15 for g = 3. For g > 4 
the bulk transition is first order. j— , 

Fisher's prediction for the bulk exponent isO 



Xb 



3-\/5 



0.19098. 



(7) 



Figure S shows a log-log plot of bulk magnetization 
vs. system size. The data correspond to different g and 
r and are all roughly parallel with slopes in the range 
0.183 — 0.187 indicated in the figure and obtained by a 
simple linear fit. In the scale of the figure oscillations are 
not visible and the estimated slopes correspond to the 
average of the approximants of Table |ll[ 



TABLE II. Approximants of the bulk exponent for various q and r. 





q 


= 2 




g = 3 




q 


= 4 


q = 6 


L 


r = 2 


r = 8 


r = 2 


r — 4 


r = 8 


r = 5 


r = 10 


r = 1.5 


5 


0.15352463 


0.18540035 


0.17119099 


0.18598351 


0.19697064 


0.19644591 


0.20327857 


0.18415244 


7 


0.15902213 


0.17126232 


0.17606416 


0.17698031 


0.17098093 


0.17386925 


0.16655933 


0.18858150 


9 


0.16285921 


0.18653590 


0.17971019 


0.18728278 


0.19301242 


0.19257056 


0.19639122 


0.18942360 


11 


0.16572102 


0.17921449 


0.18137938 


0.18203550 


0.19536097 


0.18054979 


0.17782952 


0.18957220 


13 


0.16802871 


0.18804307 


0.18299470 


0.18816428 


0.18974244 


0.19214731 


0.19569786 


0.18635074 


15 


0.16991419 


0.18254668 


- 


0.18411149 


0.19515526 


0.18319909 


- 


- 
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FIG. 3. Power-law behavior of the bulk magnetization with 
the lattice size for several values of q and r (inset shows dif- 
ferent disorder amplitudes for q = 4). The values indicated 
correspond to a simple power-law fit. 



We mention that, surprisingly, we observed a weak de- 
pendence of the effective critical exponents with the am- 
plitude of disorder r in the range [2; 20], while impor- 
tant cross-over effects with tha-atrength of the disorder 
were observed, e.g. by Picco,t3 in the case of a two- 
dimensional Potts model with a homogeneous disorder. 
This is an evidence of the strongly attractive character 
of the random fixed point. 



B. Transverse magnetization profiles 

Conformal invariance techniques are extremely ac- 
curate for studying second order phase transitions in 
isotropic pure systems. Conformal symmetry requires 
translation, rotation and scale invariance; obviously a 
single disorder realization does not have such properties, 
but it is plausible that conformal symmetry is restored 
after averaging over different disorder configurations. In- 
deed in recent investigations, conformal mappings of cor- 
relation functions and order parameter profiles have been 
successfully applied ijj-two-dimensional systems with ho- 
mogeneous disorder .t^ The case of correlated disorder, 
studied in this paper, is quite different as the system, 
even after disorder average, exhibits infinite anisotropy 
at the random fixed point (see Sec. M). Therefore con- 
formal invariance should not hold. In spite of these re- 
strictions, the magnetization profiles of the RTIC with 
fixed-free boundary conditions have been quite satisfac- 
torily fitted with the confamial expression,cj, obtained 
e.g. by Burkhardt and XueEJ 



[mi(/)]av = ^(7r/L)-^[sin7rC]-"1cos7rC/2]-S 



(8) 



where C = (' ~ l/2)/i- The success of such an ap- 
proach probably results from the geometrical constraints 
induced by the strip shape of the system combined to 
scaling arguments: In this sense, if Eq. (|8|) turns out to 
be valid also for the disordered Potts chain, one should 
conclude that such a form of profile is not a strict con- 
sequence of conformal invarianiie, which does not hold 
in anisotropic critical systems]^ but follows from more 
general grounds. 
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FIG. 4. Magnetization profiles for several values of q and 
r (inset for q = A and several disorder amplitudes) and the 
conformal expression with the Fisher exponents (solid line) 
with only the amplitude as a free parameter. 

In this Section, we indeed ask the question whether 
Eq. (^) correctly describes the transverse magnetization 



profiles in the random Potts chain. Senthil and Majum- 
darEj indeed argued that the scaling functions should be 
independent of q. The quantity ra{l)L^'' being such a 
scaling function, if Eq. (g) holds in the RTIC chain when 
L -^- CX3, it should also be valid in the same limit in the 
Potts case. Unfortunately, in our case, the largest strip 
being L = 16, strong lattice effects make the fit to ex- 
pression (g) impossible in practice (Fig. p. The same 
observation was reported for the RTIC, but these finite- 
size effects vanish for the largest systems that can be 
reached using the free fermion mapping in the Ising case. 
If Eq. (g) is valid, the quantity 



xi{L) 



ln[mL(3L/4)]av-ln[mL(L/4)], 
ln(V2-l) 



(9) 



should converge towards the surface exponent 1/2. The 
results are shown in Table III. This method cannot be 



considered as an exact determination of xi but it pro- 
vides an indication of the possible validity of (||) for this 
system. The bulk exponent xi, also follows from a three- 
point expression analogous to Eq. (0) and should be ex- 
tracted from the magnetization profile, but our data are 
not accurate enough to allow this calculation. 

TABLE III. Approximants of the surface exponent for var- 
ious q and r calculated with (Bl). The sizes in bold face cor- 
respond to chain for which there is a lattice point I such that 
{I — l/2)/L = 1/4, while the other data were obtained after 
fitting the finite-size profile. The results are compatible with 
the expected value of 0.5, seemingly supporting the expression 
of the transverse profile. 



L 


g=2-r =8 


g=3-r =4 


g = 4 — r = 5 


q=6-r=A 


4 


0.5161 


0.5195 


0.5505 


0.5709 


6 


0.5094 


0.5237 


0.5494 


0.5711 


8 


0.5039 


0.5224 


0.5411 


0.5607 


10 


0.4996 


0.5210 


0.5341 


0.5511 


12 


0.4995 


0.5202 


0.5314 


0.5422 


14 


0.5010 


0.5206 


0.5310 


0.5262 


16 


0.5000 


0.5187 


0.5269 


- 



IV. MULTI-FRACTALITY 

A. Probability distribution 

According to Fisher, the probability distribution of the 
surface magnetization to|^ for a strip of width L is ex- 
pected to behave asymptotically as: 






(10) 



This expression was successfully used for the RTIC by 
Igloi and Rieger to rescale the probability distribution of 
the surface magnetization for several strips of different 



widths L in Ref . 32 , The scahng function p was found to 



have a power-law decay in the RTIC. Forgetting about 
the oscillations, the average behavior in Fig. jsjis compat- 
ible with this observation. j— , 

The results of Senthil and MajumdaiO state that the 
scaling function p is independent of q in the thermody- 
namic limit. Here again, the widths of our strips are 
too small to give a smooth probability distribution, as 
shown in Fig. B[ The Figure shows that for the events 
which determine the average critical properties, i.e. those 
corresponding to the dominant contributions to the mag- 
netization, the probability distribution is weakly depen- 
dent on g, in accordance with universality of Eq. (10). 



On the other hand, for the rare samples for which to£ or 
m\ < 10~^, Eq (|lO)breaks down. It is however possible 
to test expression ( Iw by analyzing the moments of the 
distribution. 
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FIG. 5. Probability distribution of the surface and bulk 
(inset) magnetizations defined with 64 bins (q = 3 (squares) 
and q — 4 (circles), L — 16). 



B. Moments analysis 

In this section we analyze moments of bulk and sur- 
face magnetization distributions. This analysis will help 
to avoid the problems related to the oscillations of the fi- 
nite size estimates of the bulk exponent which, as shown 
in the previous section, hindered a reliable extrapolation 
of such quantity. Moments analysis was roeefitly used 
in studies of self-organized critical systems.l230 It was 
shown that such method provides an accurate means of 
investigation of subtle aspects that seem to emerge for 
such models. 

In our case we are interested in the scaling of the 
asymptotic behavior of the moments. For instance, in 
the case of the surface or bulk magnetization we expect: 

[{niLrU - L-^^'^K (11) 

The exponent x{^) generalizes the definitions of Eq.(ra), 



which correspond to uo — \. If the system is self- 
averaging, wp-jexpect a linear dependence of x(a;) with 
respect to w. 123 At a fixed size L, we consider the finite- 
size approximants of x('-^) defined as usual by: 



Xl{uj) = - 



ln[(m^+i)"]av~ln[(m,._i)"]3 
ln(L + 1) - ln(i - 1) 



(12) 



Plots of both quantities for L = 5, 7, ... 15 and g = 3, 
r — 4: are shown in Fig. H. The values aX lu = 1 are those 
shown in the tables | and y. The fact that Xl(0) — 
is a consequence of the normalization of the probability 
distribution. 

We notice that in both the surface and bulk case xl (^) 
for large L seems to be a constant independent of lu; this 
constant equals the surface or bulk exponent. This can 
be understood if one assumes the probability distribution 
given in Eq. (^. As a consequence of this one indeed 
obtains: 



[(to 



l) ]av 



1 

7l 
1 



himi 



(1 



dmi {mir P , /- 
V v-t^ 

1 

— , dx p 

)VlJo V(i 



In: 



^)Vl 



(13) 



where we have used the change of variables x — (to.^)'^"'"^. 
We notice from this formula that w does not modify the 
exponent of L, but enters only as a prefactor of \/L and 
consequently all the moments should have the same scal- 
ing dimension x{^)- This is what observed in Fig. y 
where the calculated xl{^) does not seem to depend on 
uj, therefore our numerical results are consistent with a 
probability distribution of the form of Eq. (p^. It is also 
important to stress that results would have been very dif- 
ferent for a probability distribution with another scaling 
variable, e.g. ml/VL- We note that in Eq. (|l^), the last 
integral seems to be nothing but a normalization, but 
as we mentioned in the previous Section the expression 
used for the probability distribution is only asymptoti- 
cally valid for the dominant events, and thus the true 
normalization of the distribution differs from Eq. (|l3|). 



?^ 




FIG. 6. Approximants of the bulk and surface magnetiza- 
tion critical exponents vs the order of the moments used to 
calculate them. The different symbols correspond to different 
sizes from 5 to 15 (g = 3, r = 4). 

An interesting consequence of these observations is 
that one could calculate surface or bulk exponents from 
an arbitrary value of lo. From the point of view of the 
finite size scaling analysis we consider values of lo where 
the dependence on the system size L is the smallest. 

For the data referring to the bulk exponent and shown 
in the figure we see that around w = 2.5 (this point is 
indicated by a thick arrow in the figure) finite size ef- 
fects are very small, all xl('^) seem to intersect into a 
single value; for which we estimate Xh — 0.1905(5) in ex- 
cellent agreement with the magnetic exponent obtained 
by Fisher (3 - ^5)74 « 0.19098 ... and also indicated by 
a thick dashed line in the figure. 

TABLE IV. Calculated bulk exponent from generalized 
scaling of the moments. 



Q = 2, r = 2 


Q = 2, r = 8 


Q = 3, r- = 4 


= 6, r = 1.5 


0.189(3) 


0.196(3) 


0.1910(5) 


0.190(2) 



Table I^ collects the values of bulk exponents calcu- 
lated with this method. For all cases studied we find a 
value of uj for which xl (w) depends only very weakly on 
L; if the intersection of the xl(ci;) for different L is sharp 
(as in the example shown in Fig. |6|) one arrives to a very 
good estimate of the exponent. 

We should also stress that, as the system size enters 
in Eq. (|l^) together with w in the form of a scaling vari- 
able (1-1- uj)y/L. This suggests that higher moments of 
the magnetizations should show weaker finite size cor- 
rections, as increasing w has the same effect of enlar ging 
the system size. This is only partially true as Eq. ( [l^ ) 
describes only the leading scaling behavior in 1/L, and 
it does not include further corrections to scaling terms. 



V. CORRELATION FUNCTIONS 

In this Section, we deal with the dynamical properties 
of the random Potts chain. We thus focus here on cor- 
relation functions in the temporal direction, which are 
defined, in the quantum formalism, by: 
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where |0) denotes the ground state. Fisher's renormaliza- 
tion group does not provide any information about this 
quantity. However, Igloi and Rieger by means of scaling 
arguments derived the following asymptotic behavior: 



This asymptotic decay of the correlation functions is un- 
usually slow with respect to usual power-law behavior 
expected at criticality. This can be understood more 
easily in the classical version of the model and it is due 
to the fact that disorder along the transfer direction is 
strongly correlatedE3 (layered disorder) and consequently 
the model displays anisotropic critical behavior. It means 
that at the critical point, distances measured in space and 
time directions are connected through power laws involv- 
ing the dynamical exponent z: i ~ P. The homogeneity 
assumption for the correlation functions thus becomes 
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where & is a rescaling factor. The choice b — t^^^ 
leads to the usual power-law decay 
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thus 
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while in the case of infinite anisotropy (z —f oo) one 
should set b = \nt and Eq. ( p^ ) follows. 

In the classical formulation of the problem followed in 
this paper, the dynamical correlation function is defined 
by: 
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where Aq and \vq) are the largest eigenvalue and the cor- 
responding eigenvector of the transfer matrix T. For the 
calculation, we used free boundary conditions on both 
edges (the magnetization defined by Eq. (Q) is thus zero 
in the whole strip) and we restricted ourselves to the spin 
in the middle of the system. Since large strip widths are 
needed in order to avoid finite size effects, average are 
now performed with ~ 25000 different samples. 
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FIG. 7. Dynamical spin-spin correlation functions for the 
RTIC {q = 2) at r = 4 for several strip widths. The number 
of samples is 25000 for L = 18 and 24000 for L = 32 and 48. 



Figure shows the correlation function calculated by 
DMRG techniques following Eq. (|l^) for g = 2 and r = 4. 
We plot the logarithm of the correlation function vs. the 
double logarithm of the distance along the time direc- 
tion (called n in Fig. nl); the dashed line is the asymp- 
totic slope predicted by Eq. (|l5|) , where we have used 
Fisher's exponent given in Eq. (ItI). This asymptotic be- 
havior seems to be in agreement with our numerical data 
for the largest system size analyzed (L = 48) for which 
we took m = 16 and 24000 disorder realizations. 
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FIG. 8. Dynamical spin-spin correlation functions for the 
g = 3 Potts chain at r = 4 for several strip widths. The 
number of samples is 30000 for L = 18 and 26000 for L = 40. 

Figure || shows the analogous plot of the time correla- 
tion function for q — 3 and r = 4. Again for the largest 
system size investigated (L = 40) one notices good agree- 
ment with the expected asymptotic limit shown as a thick 
dashed line. Differently from Fig. ^ we notice here a 
crossing of the curves representing correlation functions 
for L = 18 and L — 40. This is we believe an effect of 
the DMRG approach which does not correctly describe 
the short distance correlation functions, since it does not 
reproduce the full transfer matrix spectrum. In the other 
hand, when the number of transfer matrix products be- 
comes large in Eq. ( |l8| ) the correlation function is deter- 
mined by the extremal part of the spectrum of T, which 
is believed to be well-reproduced by DMRG methods. 
Therefore we believe that the crossing of Fig. g is due to 
the limits of the DMRG approximation, while results at 
larger n should be more accurate. 

In conclusion the analysis of the time correlation func- 
tion for (7 = 2 and q — 3 also indicates that the universal 
dynamical properties of the disordered quantum Potts 
chain are independent of q. 



VI. CONCLUSIONS 

We have provided numerical evidences supporting the 
conclusions of Scnthil and Majumdarli3 concerning the 
existence of a unique universality class for all finite value 
of the g-state Potts model perturbed by a correlated dis- 
order. For this purpose, we have measured the bulk and 
surface critical exponents and found values fully compat- 
ible with those exactly calculated by Fisher for the RTIC. 
We point out that a careful numerical study of the critical 
properties of disordered quantum chains requires an accu- 
rate average over randomness, since the most important 
source of error is due to the disorder average. This effect 
becomes essential in some circumstances, possibly lead- 
ing to wrong exponents, and may even dominate other 
contributions like finite-size effects. For that reason, we 
generated all the disorder realizations for chains of rel- 
atively small lengths, using the fact that the finite-size 
corrections turned out to be quite small (even compared 
to their importance in the corresponding pure models). 
Furthermore, we have checked the expression of the prob- 
ability distribution by an analysis of its moments: Since 
the disorder average was exact up to the accuracy of the 
DMRG method, this computation was allowed even for 
large moments. The success of this approach is largely 
due the occurrence of an optimal moment for which the 
finite-size corrections appear to be extremely weak. The 
method looks promising for further applications to dis- 
ordered systems. Finally, we have also compared the ex- 
pression of the dynamical correlation function with that 
of the RTIC. 

Another interesting feature ig the illustration, after the 
work of Aizenman and Wehr,Ll that correlated disorder 
induces a second-order phase transition in the g— state 
Potts model with q > 4. The universality class is how- 
ever very robust, in contradistinction with the homoge- 
neous two-dimensional disordered fixed point where the 
magnetic exponents continuously vary with q. 
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